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Abstract. We investigate the prospects for indirect detection of right-handed sneutrino dark matter 
at the IceCube neutrino telescope in a U{\)b-l extension of the MSSM. The capture and annihila- 
tion of sneutrinos inside the Sun reach equilibrium, and the flux of produced neutrinos is governed 
by the sneutrino-proton elastic scattering cross section, which has an upper bound of 8 x 10~ 9 pb 
from the Z' mass limits in the B — L model. Despite the absence of any spin-dependent contribution, 
the muon event rates predicted by this model can be detected at IceCube since sneutrinos mainly 
annihilate into leptonic final states by virtue of the fermionZ? — L charges. These subsequently decay 
to neutrinos with 100% efficiency. The Earth muon event rates are too small to be detected for the 
standard halo model irrespective of an enhanced sneutrino annihilation cross section that can ex- 
plain the recent PAMELA data. For modified velocity distributions, the Earth muon events increase 
substantially and can be greater than the IceCube detection threshold of 1 2 events km~ 2 yr~ 1 . How- 
ever, this only leads to a mild increase of about 30% for the Sun muon events. The number of muon 
events from the Sun can be as large as roughly 100 events km~ 2 yr~' for this model. 
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INTRODUCTION 

Supersymmetry is a front-runner candidate to address the hierarchy problem of the 
standard model (SM), and it has a natural dark matter candidate, namely the lightest 
supersymmetric particle (LSP), which can have the correct thermal relic abundance. A 
minimal extension of the SM gauge group, motivated by the nonzero neutrino masses, 
includes a gauged U{\)b-l gauge symmetry 01] (B and L are baryon and lepton number 
respectively). Anomaly cancellation then implies the existence of three right-handed 
(RH) neutrinos and allows us to write the Dirac and Majorana mass terms for the 
neutrinos to explain the light neutrino masses and mixings. 

The B — L extended MSSM alsoprovides new dark matter candidates: the lightest 
neutralino in the B — L sector [0, Hfl and the lightest RH sneutrino |@]. In this work 
we will focus on the sneutrino as the dark matter candidate^. The candidate is made 
stable by invoking a discrete 7?-parity, but in the context of a B — L symmetry, a discrete 
matter parity can arise once the U (1) b-l is spontaneously broken [6] . The B — L gauge 
interactions can yield the correct relic abundance of sneutrinos if the U(\)b-l is broken 



1 It is also possible to have successful inflation in the context of the U{\)b-l model \S\. In this case the 
dark matter candidate (the RH sneutrino) can become a part of the inflaton field and thereby gives rise to 
a unified picture of dark matter, inflation and the origin of neutrino masses 0]. 



around the TeV scale. Recently, it has been shown that it is possible to explain the 
positron excess observed in the PAMELA data OD in the context of a low scale B — L 
extension of the MSSM [0, S, OH. 

The RH sneutrino of this B — L extended model can be detected when it elastically 
scatters off a nucleus. The sneutrino-proton scattering cross section is in the 10~ H — 
10~ 8 pb range from the Z' mass limits EL ELI, large enough to be probed in the ongoing 
and upcoming dark matter direct detection experiments 0,0]. In addition, annihilation 
of sneutrinos at the present time produces LH neutrinos. It is interesting to investigate 
the possibility of indirect detection of sneutrino dark matter by using final state neutrinos 
in the IceCube neutrino telescope. This ongoing experiment plans to probe the neutrino 
flux arising from the annihilation of gravitationally trapped dark matter particles in the 
Sun and the Earth. 

THE U(\) B -l MODEL AND NEUTRINO PRODUCTION 

Since this B — L is a local gauge symmetry, we have a new gauge boson Z 1 (and its su- 
persymmetric partner). In the minimal model, we also have two new Higgs fields H[ and 
H' 2 (that are SM singlets) and their supersymmetric partners. The vacuum expectation 
values (VEVs) of these Higgs fields break the B — L symmetry [9D and result in a mass 
niz' for the Z' gauge boson. Various B — L charge assignments are allowed by anomaly 
cancelation. We choose the charge assignment in [|9D, where H' 2 couples to the RH neu- 
trinos and gives rise to a Majorana mass upon spontaneous breakdown of the U{\)b-l- 
Choosing these Majorana masses in the 100 GeV — 1 TeV range, we have three (dom- 
inantly RH) heavy neutrinos and three (dominantly LH) light neutrinos. The masses of 
the light neutrinos are obtained via the see-saw mechanism. 

A natural dark matter candidate in this model is the lightest sneutrino N. We consider 
the following two cases: 

• Case 1: A generic case where a solution to the positron excess observed by 
PAMELA is not sought. In this case the dominant annihilation channels are the 
S- wave processes AW — > AW and N*N* — > N*N* via f -channel exchange of Z'. 
There are also AW* — > AW* , // annihilation modes via ^-channel exchange of a Z' 
ox B — L Higgs fields, but these are P-wave suppressed and can be completely ne- 
glected (particularly at the present time). In this case the annihilation cross-section 
has the nominal value ~3x 10~ 26 cm 3 /sec (dictated by thermal freeze out) at all 
times. The RH neutrinos produced from dark matter annihilation quickly decay to 
LH neutrinos and the MSSM Higgs. Assuming that the mass difference between 
the RH sneutrinos and RH neutrinos is small, the RH neutrinos are produced non- 
relativistically, and hence each LH neutrino and Higgs receives an energy equal to 
half of the sneutrino mass. This produces a delta function in the energy of the LH 
neutrinos at one-half the mass of the sneutrino dark matter. A small portion of muon 
neutrinos from this initial annihilation state are scattered via neutral current interac- 
tions inside the Sun to lower energies. This produces a slight bump in the neutrino 
spectrum at low energies. As expected for monochromatic neutrinos, the spectrum 
of muons from charge current interactions in case 1 has a linear dependence on 



energy. 

• Case 2: In this case the PAMELA puzzle is addressed via Sommerfeld enhance- 
ment of sneutrino annihilation at the present time [8D. The dominant annihilation 
channel is N*N — > 00 via the s-channel exchange of the new scalar Higgs fields, 
the t or w-channel exchange of a N, and the contact term |/V| 2 2 . The cross section 
for annihilation to the 00 final state at the present time is required to be 3 x 10~ 23 
cm 3 /sec in order to explain the PAMELA data. Sufficient Sommerfeld enhance- 
ment is obtained as a result of the attractive force between sneutrinos due to the 
exchange provided that the mass of is small (< 20 GeV). The subsequently 
decays into fermion-antifermion pairs very quickly via a one-loop diagram, and it 
mostly produces t + t final states by virtue of the fermion B — L charges [8]. RH 
neutrinos constitute about 10% of the annihilation final states, while 00 compose 
the remaining 90% of the branching fraction. This branching fraction is neces- 
sary to provide a high enough leptonic particle rate to fit the PAMELA data. For 
4GeV < < 20 GeV, the final states are mostly taus (74%) and b quarks (16%), 
where the dominance of tau final states is a result of the fermion B — L charges. The 
LH neutrinos in this case arise from the three-body decay of taus and bottom quarks 
and are spread in energy signal. The delta function from the neutrino channel at the 
detector is subdominant to the other annihilation channels. The spectrum of muons 
is softened compared to case 1. 

Both the case 1 and case 2 scenarios of our model display a crucial signature difference 
when compared to the standard neutralino LSP in the MSSM. The energy distribution 
of the produced LH neutrinos from the RH neutrino decay is a delta function occurring 
at half of the sneutrino mass. Other annihilation channels in this model, as well as those 
available in the MSSM, produce additional neutrino signal via three-body decays such as 
T~ — > e~ v T v e . This difference opens up a significant possibility to differentiate between 
the B — L model and the MSSM with the help of the differential energy spectrum of the 
detector event rates. 



ICECUBE NEUTRINO SIGNAL RESULTS 

Sneutrino annihilation in the Sun and the Earth produces an expected neutrino flux and 
charge current interactions in IceCube. The muon track signal can be differentiated from 
the cosmic ray induced background by selecting upward-going and contained muon 
events, subtracting atmospheric neutrino background expectations and making angle 
cuts on the muon tracks. The muon and neutrino fluxes are modeled by calculating 
the number of gravitationally captured sneutrinos and then considering the propagation 
and detection of the produced neutrinos. Using DarkSUSY and WimpSim tables, our 
calculations account for neutrino oscillation, loss via charged current interactions and 
scattering via neutral current interactions [15, DarkSUSY default parameters are 



used. For both case 1 and case 2, the maximum spin-independent cross section allowed 
by the Z' limits is used. We assume equal branching to the three flavors of LH neutrinos, 
since the results of [9] do not depend critically on the choice of neutrino flavor branching 
ratios in either case. 
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FIGURE 1. Total muon rates detected at the Earth from annihilation of sneutrino dark matter in the Sun 
as a function of the sneutrino mass. The results are for one year of detection with IceCube. Case 1 (case 
2) is the highest (lowest) peaked line. The dotted line denotes the mass range where one cannot explain 
the PAMELA data using case 2 anymore. 

Fig. Q] shows our results for the total muon rate integrated over energy as a function of 
the sneutrino mass for annihilation in the Sur0. 

The general decrease of the event rates for higher is reflective of the decrease 
of the neutrino flux due to the kinematic suppression of sneutrino capture (~ l/m^ for 
large masses). The linear increase at low is explained by the linear dependence of the 
cross section for charged current interactions on the energy of neutrinos at the detector 
(~ m^). LH neutrinos are produced in two-body decays in case 1 versus three-body 
decays in case 2, and hence have a higher energy. Thus at lower values of sneutrino 
mass, the cross section for conversion of neutrinos to muons at the detector is larger in 
case 1. However, for large sneutrino masses case 1 has a smaller signal than case 2 since 
neutrinos get absorbed via charged current interactions or lose energy via neutral current 
interactions inside the Sun more efficiently because of their larger energy. 

According to refs. Il3l [Till, in the case of the Earth more than 12 events are needed 
for a DM mass between 70 GeV and 4 TeV. In the case of the Sun the number of events 
needed drops linearly as a function of mass starting from 300 events at 70 GeV down 
to 70 events at 300 GeV. Beyond 300 GeV up to 4 TeV, the number of events needed 
remains fixed at 70. This provides a hint that one could detect the event rates caused by 
sneutrinos despite some differences between the sneutrino and neutralino dark matter 
spectra used to calculate these sensitivities. Hence, it might be possible to detect muon 
neutrinos produced by sneutrino annihilation for sneutrino masses around 300 GeV for 
the Sun, cf. Fig. [IJ Note that a large range of masses would be accessible with only an 
order of magnitude improvement in sensitivity. 

For a known sneutrino mass, observation of a muon signal exceeding the number 
given in Fig. Q] will rule out the B — L model. The largest number of muon events from 
the Sun in the entire depicted mass range is 58 km~ 2 yr _1 (36 km~ 2 yr _1 ) for case 1 
(case 2). 



The apparent discrete nature of these plots occurs because only a few values of sneutrino mass are 
recorded in the WimpSim tables used by DarkSUSY; the program interpolates between these points. The 
effect is numerical and not physical. 



DARK MATTER DISC IN THE MILKY WAY 



In our analysis, we assumed a Gaussian-like velocity distribution for dark matter parti- 
cles with a three dimensional velocity dispersion of o v = 270 km sec -1 and \vs U n\ = 220 
km sec -1 for the velocity of the solar system with respect to the halo. However, there 
are recent speculations about the existence of a dark matter thick disc in the Milky Way 
in addition to the baryonic one, see e. g. |17|, HI]. This dark matter disc would have dif- 
ferent ranges for the solar system velocity and velocity dispersion: \vs m \ ~ — 150 km 
sec -1 and o v ~ 87 — 156 km sec -1 lfl9h . A modified velocity distribution substantially 
enhances the Earth muon rate for the sneutrino dark matter beyond the detection thresh- 
old of 12 km~ 2 yr _1 over a large part of the velocity parameter space [9]. It also raises 
the maximum Sun muon rate to 78 events km~ 2 yr~ ! (48 km -2 yr _1 ) in case 1 (case 2). 



COMPARISON WITH MSUGRA 

Minimal supergravity (mSUGRA) is a constrained version of the MSSM with the light- 
est neutralino as its dark matter candidate. In this section we compare the B — L model 
with the mSUGRA hyperbolic branch/focus point mSUGRA scenarios, where the dark 
matter has a large Higgsino component and mo is very large. The neutralino has a large 
capture rate in this region due to a large Higgsino component that results in a large 
spin-dependent scattering cross section via Z exchange. 
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FIGURE 2. Total Sun-annihilation muon rates inside the detector for mSUGRA hyperbolic 
branch/focus point scenarios as a function of the neutralino mass. The results are for one year of detection 
with IceCube. 

Fig.|2]shows the total Sun muon rate as a function of the neutralino mass for mSUGRA 
hyperbolic branch/focus points. A comparison with Fig. Q] shows that these scenarios 
always have a higher total muon rate in the plotted mass range than the B — L model. 
The hyperbolic branch/focus point models yield larger muon rates by between more than 
one order of magnitude and a factor of 1.5 for dark matter masses in the 100 — 800 GeV 
range. These higher rates are explained by the spin-dependent scattering cross sections, 
which are a few orders of magnitude larger than the limits on the spin-independent cross 
section for the sneutrino dark matter. The spin-dependent scattering cross section for 
the B — L model is zero because U(l) b-l is a vectorial symmetry. Since the Sun mainly 
consists of hydrogen, the spin-dependent piece contributes dominantly for the mSUGRA 
case. 




However, it is interesting that despite having a much smaller scattering cross section, 
the B — L model can yield muon rates that are roughly comparable to the mSUGRA 
scenarios. Sneutrino annihilation dominantly produces leptons, i.e., RH neutrinos in 
case 1 and taus in case 2, which subsequently decay to LH neutrinos 100%. On the other 
hand, neutralino annihilation in the hyperbolic branch/focus point scenarios dominantly 
produces quark final states that have a small branching ratio for decay to neutrinos. 
Furthermore, despite lower event rates, sneutrino dark matter still produces a distinctive 
linear spectrum in the muon flux. This feature is caused by the delta function in energy 
for the neutrino spectrum and can be used to distinguish between the B — L model and 
the hyperbolic branch/focus point scenarios as long as energy binning of the differential 
muon rate with respect to the energy is precise enough at IceCube. 
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